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Abstract 

The tolerance of bacterial populations to biocidal or antibiotic treatment has been well documented 

in both biofilm and planktonic settings. However, there is still very little known about the mechanisms 

that produce this tolerance. Evidence that small, non-mutant subpopulations of bacteria are not affected 

by antibiotic challenge has been accumulating and provides an attractive explanation for the failure of 

typical dosing protocols. Although a dosing challenge can kill all the susceptible bacteria, the remaining 

p^ persister cells can serve as a source of population regrowth. We give a robust condition for the failure of 

(— ^ a periodic dosing protocol for a general chemostat model, which supports the mathematical conclusions 

^\j and simulations of an earlier, more specialized batch model. Our condition implies that the treatment 

I protocol fails globally, in the sense that a mixed bacterial population will ultimately persist above a 

^ level that is independent of the initial composition of the population. We also give a sufficient condition 

'""5 for treatment success, at least for initial population compositions near the steady state of interest, 

^SJ corresponding to bacterial washout. Finally, we investigate how the speed at which the bacteria are 

1-H wiped out depends on the duration of administration of the antibiotic. We find that this dependence is 

not necessarily monotone, implying that optimal dosing does not necessarily correspond to continuous 

administration of the antibiotic. Thus, genuine periodic protocols can be more advantageous in treating 

a wide variety of bacterial infections. 
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1 Introduction 



The failure of antibiotic treatments to eliminate bacterial infections has become both more evident and 
^~* better understood in the past several decades. Although there is evidence that the over use of antibiotics 

Jzl has amplified the number of chromosomal-resistant bacteria [18' , it is becoming increasingly clear that there 

^b are other mechanisms that protect populations of bacteria. Many of these mechanisms depend on whether 

Q^ the bacteria exist in a biofilm or not [H [101 [HJ H]. In particular, the notion that small sub-populations 

^~H of bacteria may display innate tolerance to various biocides has been proposed as a possible reason for the 

f~^ failure of treatment for bacterial infections [16, 2, 14, 5 . Bacteria within a biofilm are enmeshed in a physical 

f^ gel that provides a secondary boundary that may allow small numbers of bacteria to evade the antibiotic; 

OO therefore, the failure to eliminate the entire population can allow the population to regrow. 

It should be noted that populations of planktonic bacteria also contain these highly tolerant or persister 
J> cells [9l [21] . Thus understanding the process of persister formation and the response of the population to 

biocidal application is fundamental to developing dosing protocols and treatments in both batch culture and 
biofilm populations. 
C^ As in many areas of biology, mathematical modeling has been used as a counterpart to experimental 

observations. Because there are several hypotheses regarding the mechanism of persister formation, math- 
ematical modeling can be used to provide insight into the success of failure of treatment protocols as well 
as the consistency of various hypotheses. Currently, there are at least two distinct hypotheses concerning 
persister formation-one of which is senescence. In this case, persister cells are assumed to be those that have 
undergone many division cycles. It is known that asymmetric division leads to degradation of parts of the 
cellular machinery that may be the underlying cause of persistence [12] • Mathematical analysis of a model 
of senescence has been described in both chemostat and biofilm settings [TH [T] . 

A different hypothesis argues that persisters are a phenotype that is expressed at a rate that depends 
on the external environment. Although the biological details are not well understood, it is thought that 
this might be due to toxin/ antitoxin interaction or some other stress response [Tni [131 [H [5TJ [T7]. This has 
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been investigated mathematically as well [21 [3 [6l [12] . In [5], a very simple model of persister formation was 
developed and optimal dosing protocols, that entail alternating application and resting, were described. In 
[B], toxin/antitoxin interaction was explicitly included and the resulting model was analyzed in a chemostat. 
In both of these investigations only one particular form of the kinetics was analyzed. One of the goals of 
the current investigation is to extend these results to a more general form for kinetics. This is an important 
process, because it has been shown that not only is there a successful dosing protocol, there is an optimal 
protocol. Without knowing how robust the model is, there is no way to determine how robust the conclusions 
are. Here we take a model of the form proposed in [S] and give a condition that yields successful dosing 
protocols, at least for some initial makeups of the population. 

The manuscript is organized as follows: We begin by describing the model for the dynamics of the bacterial 
population in response to antibiotic challenge. We then develop the theory by analyzing two extreme cases 
(no dosing and constant dosing) and the intermediate case. This leads to a local sufficient condition for 
treatment success. Next we give a condition for global treatment failure, supported by numerical simulations 
of the model. We also show that the speed of eradication does not necessarily depend monotonically on the 
duration of the administration of the antibiotic. 

2 Model 

Consider the following chemostat model: 

Bs - [{l-kd{t)~ki{t))f{S)-D]B, + kg{t)Bp (1) 

Rp = ki{t)f{S)Bs-[kg{t)+D]B.p (2) 

S = DiS°-S)-^^^^ (3) 

where Bs is the concentration of the cells which are susceptible to antibiotics, Bp is the concentration of 
the persister cells which are not affected by the antibiotic and S is the concentration of the nutrient. This 
model deviates from the one in [5] because it is a chemostat model, which is reflected in the additional loss 
terms at rate D (called the dilution rate or washout rate) , and the inflow (at the same rate D) of nutrient 
with an input concentration S^. The per capita growth rate of the susceptible cells is denoted by f{S), for 
which we assume the following throughout the rest of the paper: 

/ : M+ — > M+ is smooth and increasing and /(O) = 0. 

The persister cells do not consume nutrient, hence the lack of a corresponding growth term in the Bp- 
equation. The conversion of nutrient into new biomass occurs with a yield oi Y Cz (0, 1). 

The remaining functions kd{t),ki{t) and kg{t) are non-negative, time-varying functions which describe 
the effect of antibiotics on the population. First, kd{t)f{S) is the killing rate of the susceptible population. 
Note in particular that the killing rate is proportional to the growth rate of the cells. It is positive when 
both antibiotic and nutrient are present, but zero when either one is missing. Secondly, ki{t) is the rate at 
which susceptible cells switch to persister cells. It is also positive when antibiotic is present, but zero when 
this is not the case. Finally, kg{t) is the rate at which persister cells revert to the susceptible state when 
antibiotic is absent (and zero when antibiotic is present). 

Since antibiotics are administered to the reactor vessel in a controlled (lab) environment, we will make 
the simplifying assumption that the functions kd(t),ki{t) and kg(t) are r-periodic (for some given r > 0), 
and of the bang-bang type with simultaneous switching instances: For some p e [0,1], and for positive 
parameters kd,ki and kg, there holds that 

Mi) = |'^'°''^tO,pr) ^ ,^^,^^ikiorte[0,pr) ^ ..^ M^) = ^ ^"' ^ ^ ["'^"^ . (4) 
^ ' [0 for te[pT,T) ^ ' [0 for t G [pr, r) ' [h for t G [pr, t) ^ ' 

Thus, antibiotics are present during a fraction p of the period r, and absent during the remaining fraction 
1 — p of the period. 

Clearly, this is a simplification of reality because the concentration of an antibiotic is not expected to 
be of the bang-bang type. In a more realistic model, the functions kd{t),ki{t) and kd(t) would be replaced 
by functions depending on (at least) a new state variable for the concentration of the antibiotic, and the 



periodicity would arise through a periodic forcing term in the equation for this new variable. We leave the 
study of such a model to the future. 

Throughout the rest of this paper we assume that the net effect of the antibiotic alone is lethal to the 
susceptible population: 

l-kd-ki<0. (5) 

Note that this assumption is valid for the parameter values related to the experiments described in fS]. 

The main purpose of this paper is to investigate how the behavior of system ([l]) — ([3| with Q changes 
qualitatively in terms of p. 

3 Preliminary results 

In this section we collect a couple of basic results concerning the dynamical behavior of (IT]) — (l3| with Q . 
For a real- valued function x{t), we denote the extended real numbers liinhift^oc x(t) and limsupj^^^ a;(i) 
by Xoo and x°° respectively. 

Lemma 1. System (jll) ~ (|3| with Ml has M.^ as a forward invariant set, and it is dissipative. 

Proof. The first assertion is obvious. Dissipativity follows by considering the dynamics of 

M = Bs + Bp + YS, 

given by 

M = £1(^5'° - M) - kd{t)f{S)B, < D{YS° - M), 

and hence 

M°° < YS". 

n 

Not only are all state components of every solution ultimately bounded from above by some constant 
which does not depend on initial conditions, we also notice in the following Lemma, that S{t) is ultimately 
bounded from below by some positive constant which is independent of initial conditions. 

Lemma 2. There is a constant 9 > such that for all solutions o/ (fl]) — ([3| with (H), there holds that 

Proof. Consider the function g{x) :— f (x) S° — D{S^ — x) . Then g is increasing with ^(0) < and g{S°) > 0, 
hence by the intermediate value theorem, there is a unique 9 G (0, 5°) such that g{9) = 0. We will show 
that S'oo > 9. If not, then since B'^ < YS^ by the proof of Lemma ll] it follows from Corollary 2.4 in [23] 
-a consequence of the famous Fluctuation Lemma- applied to ^ that 



> lim inf 

t — >oo 



i^(50_5^)_/(^-)^«W 



Y 

> D{S''-Soo)-f{Soo)S'' 

> D{S° -9)- f{9)S°, 

which contradicts that g{9) = 0. D 

4 Analysis of the extreme cases p = and p = 1. 

First we study the cases where antibiotic is either present or absent for all times. Our conclusions are the 
expected ones: When antibiotic is present continuously, all susceptible cells are killed, and consequently 
the persisters die out as well since they can only survive if susceptible cells become persisters. When the 
population is never exposed to antibiotics, then all persisters revert to the susceptible state, and ultimately 
the population will consist entirely of susceptible cells (provided the dilution rate is not too high). 

Lemma 3. 1. Case p = I (continuous antibiotic dosingj. The steady state {Bg, Bp,S) = (0,0,5"°) 
of system (fl]) — ([3| with m is globally asymptotically stable. 



2. Case p — (never antibiotic dosingj. If D < f{S^), then all solutions {Bs{t),Bp{t), S{t)) with 
-Bs(O) + Bp{0) > converge to the steady state {Bs, Bp,S) — {B*,0,S*), where S* G (0,5*") is the 
unique positive value satisfying f{S*) = D, and B* = {S^ — S*)Y . 

Proof. 1. If p = 1, then kj^{t) = k^, ki{t) = ki and kg{t) = for all t. By (Is]) it follows that for all 
solutions, Bs < —DBg, and hence Bs{t) ^ as t ^ oo. This suggests that we should study the linear 
limiting system 

Bp — —DBp 

S ^ D{S'^-S) 

whose solutions clearly converge to {Bp,S) — (0, S*"). The conclusion now follows immediately by 
applying Theorem F.l in I20j . 

2. If p — 0, then fcrf(i) = fc;(t) = and kg{t) = kg for all t. Then the restriction on the initial condition 
implies that Bs{t) > for all i > 0. Notice that Bp{t) ^ as t ^ oo, suggesting we should study the 
limiting system 

B, = [f{S)-D]Bs 



Y 

This is the classical chemostat model of a single cell population growing on a single nutrient, described 
in [2U]: all solutions with i?s(0) > converge to {Bs,S) = {B*,S*) under the assumptions of the 
Lemma. The conclusion follows again by Theorem F.l in pp] . 

n 
5 The case of periodic dosing: p E (0, 1). 

In this section we deal with the r-periodic model (III — ([3]) with Q, assuming that p e (0, 1). We will also 
assume that the dilution rate is not too large, as in Lemma [3] 

D < /(5°). (6) 

Notice that system (fl} — S with Q has a steady state Eq '■— {Bs,Bp, S) — (0, 0, 6*0), regardless of the 
value of p. To determine its stability properties, we determine the r-periodic variational equation at Eq: 

^{l-kd{t)-ki{t))f{S°)^D 
f{S')ki{t) 

/(g°) 

Y 

and the Floquet multipliers of this system are e~^'^ (which is of course inside the unit circle of the complex 
plane) and the Floquet multipliers of 

{l~kd{t)~kl{t))f{S'>)~D kg{t) 

f{S°)kl{t) -{kg{t) + D) 

Using Q , these Floquet multipliers are the eigenvalues of the following matrix: 

$ _e(l-p)rA2gPTAl^ (7) 



kg{t) 


■ 


{kg{t) + D) 








-L> 



{l-kd-ki)f{S^)-D 0\ „^^ ,^._ff{S")-D kg 



where 

"^'-■^['^ '^\J{S^Y ' ^ -"z^j'^^dA^:^!^-'^-^ - Hkg + D)) (^) 

are triangular, quasimonotone matrices. Notice also that Ai is Hurwitz by ([5|, and that A2 has one negative 
and one positive eigenvalue by ^. 

Since Ai and A2 are quasi-monotone (i.e. their off-diagonal entries are all non-negative), it follows that 
their matrix exponentials are (entry-wise) non-negative, triangular matrices and then their product $ is a 
(entry-wise) positive matrix whose spectral radius p(<i>) is an eigenvalue by the Perron-Frobenius Theorem[3]. 
Consequently, to determine stability of Eq, we need to establish whether or not p{^) is inside the unit circle: 
If p(<i?) < 1, then Eq is locally asymptotically stable. If p{^) > 1, then Eo is unstable. Summarizing, we 
have estabhshed 



Theorem 1. Let p G (0, 1), and assume that im holds. Then the steady state Eq = (0, 0, S'^) is locally stable 
for ^-^ with ^ ifp{<i>) < I, but unstable if p{<^) > 1. 

Our main concern is knowing how p(<&) varies as a continuous function of p (this variation is continuous 
since eigenvalues of a matrix are continuous functions of its entries, and clearly the entries of $ are continuous 
in p). For p = (never using antibiotic), and hence also for p near by continuity of p{^), we have that 
p(<I>) = p(e'^'^^) = e^^*^'^ )-D)r ^ -^ This is in accordance with Lemma pi where it was shown that all 
solutions with Bs{0) + -Bp(O) > converge to {Bl,0,S*), and thus Eq = (0,0,5°) must be unstable. For 
p = 1 (using antibiotic continuously), and hence also for p near 1, we have that yo(<&) = p(e'^^^) = e"-^'^ < 1. 
This is in accordance with Lemma l3] as well because it was shown there that all solutions converge to Eq in 
this case. 

We can actually determine the dependence of p($) on p explicitely because fortunately, both Ai and 
A2 are diagonalizable (their eigenvalues are distinct), which simplifies the computation of their matrix 
exponentials somewhat. It is easily verified that 



where 



Ai = TiDiT^^ and A2 = T.2D2T^\ 



rj._:i-Ki-ki 0\ ^ _({l-ka-ki)f{S'')-D 
^'-^ ki ij ' ^1 - I ^D 



and 



and thus using (|7| that 



^-^0 -(/(50) + fc,)y' '"I -(fcg + z?);' 



A lengthy algebraic calculation shows that the positive matrix *& = K^ » is given by 

\^<1>21 1'22/ 

cc&n = -(l-fc,-fcO(/(5°) + fcff)e^[^(^'')-^-P('=-+'=')/(^")l 

+kgki (e(i"P)''(-''('5°'~-°' -e~(i~P'''('=''+^)) (e"P^-° -eP^I^^^''''^'''^^^'^''^^-"!) 

C$21 = kiifiS") + kg) e-(i-p)-(fe«+^) (e-P-D ~ eP-[(i-fc^-fcO/(s°)-i?]^ 
C$22 = -if is') + kg){l - fc, - fc,) e--«i-P)'=«+^) 



and 
where 



C=-{l-kd-kl){f{S°)+kg) 

is a positive constant, independent of p. 

Since $ is a positive matrix, its spectral radius can now be calculated explicitly in terms of its entries: 

,^. _ ^11 + ^22 + v/(^ll ~ ^22)^ + 4^12^21 , . 

6 Conditions for treatment failure 

In this section we show that the spectral radius p{^) also plays a key role in the global behavior of system 
(jTl — (|3| with (|4| and p G (0,1). We will show that if p(<&) > 1 , then not only is Eq unstable as we have shown 
in Theorem [T] but treatment fails globally, because both cell populations persist uniformly. In addition we 
will show that there are positive periodic solutions. 

Theorem 2. Let p e (0, 1), and assume that ([6| holds. Lf p{^) > 1, then treatment fails and the population 
is uniformly persistent, i.e. there is some e* > (independent of initial conditions), such that all solutions 
0/ (IT]) — ([3| with Q and Bs{0) > 0, have the property that: 

Bs{t) > e*, and Bp{t) > e*, for all sufficiently large t. 

Moreover, there are r-periodic solutions {Bs{t), Bp(t), S{t)) with B^it), Bp(t) > for all t. 



Proof. Define the following matrix: 



|>(e) = e(l~P)'^'^2(c) gPriiCe) 



where 



M.) :- f^<*° - '> -i'^.t^'lT + '' - " _y , and AM .- (^<*" " '> " ^ 







-(kg+D) 



V kifiS^ - e) 

Notice that l>(0) = <f>, and thus since p((f>) > 1, it follows that 

p($(e)) > 1, for all sufficiently small e > 0, (10) 

as well, because the spectral radius of any matrix is continuous with respect to its entries. We fix some e > 



such that (10) holds 



We will first show that Bs is uniformly weakly persistent, i.e. that there is some e' > such that if 
Bs{0) > 0, then B"^ > e' . By contradiction, if Bg is not uniformly weakly persistent, then there is some 
solution {Bs{t),Bp{t),S{t)) with 5^(0) > such that 



BT < 



YD 



(11) 



By Corollary 2.4 in [23J apphed to ([s]), and since S°° < 5*° by ([s]), we have that 

f{Soo)Bs{ty 



> lim inf 

t^OQ 



i?(5° - S^) - 



Y 



> DiS'^-S^) 



/(^°)i?g 
Y 



and hence by (11) that 



2 

Thus, for some T* > 0, there holds that S° - e < S{t) < S° + e for all t > T* . It follows from ([l]) - (H), that 
for all t>T*: 



Bs\ . ff{S° - e) - {kd{t) + ki{t))f{S'> + e)-D 
bJ-[ ki{t)f{S^-e) 



kg{t) 

-{k,{t)+D: 



(12) 



where the vector inequalities should be interpreted componentwise. Notice that the vector field on the 



right-hand side of ( 12 1 is that of a r-periodic, cooperative linear system whose principal fundamental matrix 



solution evaluated over one period r equals <&(e). By Kamke's comparison Theorem (see e.g. Theorem B.l in 
Appendix B of [20]) it follows that for all t >T* , the vector {Bs{t), Bp{t)Y' is not smaller (component-wise) 
than the solution starting in {Bs{T*), Bp{T*))^ of the r-periodic, cooperative linear system with vector field 
given in the right-hand side of (12). But all non-zero, non-negative solutions of the linear system diverge 



because p(<i>(e)) > 1. Then so does {Bs(t), Bp(t)), and this contradicts (11 1. We have thus shown that Bs is 
uniformly weakly persistent. 

Next we establish that Bg is in fact uniformly strongly persistent. This follows from Theorem 1.3.3 
in [24], applied to the map P which maps (Bs(0),Bp(0), S'(O))'^ e X to (B,.(r),Sp(T), ^(t))^, where X := 
{{Bs, Bp, S)^ e Rl I Bs+Bp+YS < YS° }, Xo := {(B„ Bp, Sf £ X \ B^ ^ 0} and OXq := {(B„ Sp, Sf e 
X \ Bs = 0}. The map P is continuous and maps Xq into itself, and it has a global attractor because it is 
compact and dissipative. It follows that there is some e^ > 0, independent of initial conditions, such that if 
Bs{0) > 0, then liminf„^oo Bs{nT) > el, and also that liminft^oo Bs{t) > ej by Theorem 3.1.1 in [24] . 

Next we show that uniform strong persistence of Bs, implies uniform strong persistence of Bp. Consider 
equation ^ and notice that for all sufficiently large t: 

Bp > ki{t)f{e)Bs - {kg{t) + D)Bp > ki{t)f{0f^ - {kg{t) + D)Bp, 
where 9 is the positive constant from Lemma |2] It is not hard to show that the linear equation 

z = ki{t)!{e)'^-{k,{t) + D)z, 



has a positive r-periodic solution p{t) and that all non-negative solutions converge to it. Therefore, it follows 
that for all sufficiently large t, 

B,it) > ^, 

establishing uniform strong persistence for Bp, since Poo is independent of initial conditions. We conclude 
the proof of uniform strong persistence of Bg and Bp by setting e* = min{ei, ^}- 

Finally, to show that there are r-periodic solutions with Bg{t),Bp{t) > 0, we apply Theorem 1.3.6 from 
[24 applied to the continuous map P defined above. We have already remarked that this map is continuous, 
maps Xq into itself, is dissipative and compact, and we have just proved that it is uniformly strongly 
persistent with respect to {Xo,dXo). Observe also that Xq is relatively open in X, and that Xq is convex. 
Then by Theorem 1.3.6 from |24j . the map P has a fixed point in Xq, and this in turn implies the existence 
of T-periodic solutions with Bs (t) > for (IT|) — (l3| . The same argument as the one used above to establish 
uniform strong persistence of Bp(t), shows that these r-periodic solutions are such that Bp(t) > as well. 

D 

7 Numerical example 

We use the following numerical values: The per capita growth rate is of Michaelis-Menten type: 



/i5 
and the parameters 



f{S) = 7"^^, where ^ = 0.417 hs"^ and k^ = 0.2 mgr\ 



ka ^3, ki = 0.1, kg ^ 0.5 h"\ and r = 10 h. 
The chemostat setting requires that we specify two additional parameters: 

D = 0.1 h"^ and S° = 1 mgr\ 

It can be easily verified that these choices satisfy the conditions (Is]) and (|6l. Note that there is no need to 
specify the yield coefficient Y in order to calculate yo($). Indeed, the matrix <& in ( [7| d oes not depend on Y. 
The graph of the spectral radius p(<&) in terms of p, determined using formula (|9|, is given in Figure fl] 
We see that p{^) = 1 when p is approximately equal to 0.242. Clearly, p{^) is not monotone, and it has 
a global mimimum of approximately 0.098 which is achieved at p ~ 0.612 (determined numerically using 
Mathematica) . Thus, the optimal strategy in a dosing experiment with a period of r = 10 hours occurs for a 
dosing duration of pr « 6.12 hours. Here, optimality means that eradication happens as quickly as possible. 

Let us also illustrate what happens if p equals approximately 0.205. Then p{^) equals approximately 
1.436, implying that treatment fails. It appears that the solution of (IT|) — pi with Y — 1 (and all other 
parameters as above) starting from the initial condition {Bs,Bp, S) — (0.3,0,0.4) converges to a r-periodic 
solution, see Figures |2] and [3] These observations are in accordance with Theorem |2] 

We also see that ii p —^, then p{^) equals approximately 0.426. Then it follows from Theorem [l] that 
Eq is locally stable. Figures |4] and |5] illustrate this by indicating that for the solution with the same initial 
condition {Bs,Bp, S) = (0.3,0,0.4) as above, treatment is successful. 

Finally we remark that the spectral radius of <& may be monotone: If we change the value of kd from 3 to 
1, and leave all other parameters unchanged, then p{^) is a decreasing function of p G [0, 1], and achieves its 
minimum at p = 1. This indicates that for this case, the optimal strategy is to use antibiotics continuously, 
see Figure [6] 

8 Discussion 

Bacterial infections are a source of problems in a wide variety of situations including industrial, environmental 
and clinical settings. Growing understanding of the inability of antibiotics and biocides to treat these 
infections has driven investigations into the cause of the failure of treatments. It is becoming increasingly 
evident that persister cells must play an important role in protecting populations of bacteria. It has been 
observed that other protective mechanisms including physiological and physical processes are not sufficient to 
explain the observed failures [HH]. Moreover, because it is very difficult to investigate and classify persister 
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Figure 1: Spectral radius of $ for p e [0, 1] (kd = 3, other parameters in text.) 
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Figure 2: Times series for Bg and Bp. {p w 0.205, p($) ~ 1.436 so treatment fails) 
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Figure 3: Times series for S and Bs + Bp. {p « 0.205, p($) ~ 1.436 so treatment fails) 
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Figure 4: Times series for Bs and Bp. (p = |, p((f>) « 0.426, so treatment succeeds) 
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Figure 5: Times series for S and Bg + Bp. (p = ^ p{^) « 0.426, so treatment succeeds) 
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Figure 6: Spectral radius of $ for p G [0, 1] (kd = 1, other parameters in text). 



cells experimentally, mathematical modeling can play an important role in supporting hypotheses as well as 
generating useful predictions. 

We have described and analyzed a general model for the dynamics of persister formation in response to 
antibiotic challenge. We have been able to provide a condition for the success/failure of antibiotic challenge 
in a chemostat. We have found that there is an optimal strategy, that is one that kills the bacteria the 
fastest. These results indicate that periodic dosing is an effective treatment protocol for a variety of bacteria, 
substantially strengthening the results in 5 . The theoretical results were also confirmed by direct numerical 
simulations. 
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